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Abstract 

We calculate the hadronic light-by- light contributions to the muon g—2. 
We use both 1/iVc and chiral counting to organize the calculation. Then 
we calculate the leading and next-to- leading order in the l/N^ expansion 
low energy contributions using the Extended Nambu-Jona-Lasinio model 
as hadronic model. We do that to all orders in the external momenta and 
quark masses expansion. Although the hadronic light-by-light contribu- 
tions to muon g — 2 are not saturated by these low energy contributions we 
estimate them conservatively. A detailed analysis of the different hadronic 
light-by-light contributions to muon i^r — 2 is done. The dominant contri- 
bution is the twice anomalous pseudoscalar exchange diagram. The final 
resuh we get is ajjg^*-'^y-"s^* = (-9.2 ± 3.2) • 10"^°. This is between two 
and three times the expected experimental uncertainty at the forthcoming 
BNL muon g — 2 experiment. 
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1 Introduction 



The forthcoming experiment at Brookhaven National Laboratory (BNL) ||T| plans 
to measure the anomalous magnetic moment of the muon = (gf^ — 2)/2 with an 
accuracy around ± 4 • 10~^°, improving by a factor more than twenty the previous 
experimental determination at CERN 0: 

a;fP = 11659 230(84) ■ 10"^°. (1.1) 

This expected impressive performance has motivated the recent raised interest in 
obtaining a more accurate theoretical prediction of within the Standard Model 
(SM), for reviews see 0. One of the reasons is that with a theoretical uncer- 
tainty of the same order as the aimed BNL uncertainty, the anomalous magnetic 
moment of the muon could become a precision test of the quantum corrections 
of the electroweak sector of the Standard Model. The one-loop electroweak con- 
tributions to are @] 

af^ione - loop) = 19.5(0.1) ■ 10"^°. (1.2) 

These are order ni?Gp with Gp the Fermi coupling constant and m the muon 
mass. Electroweak contributions from two- loops of order m?Gp{a/TT) have 
been recently reanalyzed by [P, 0] for the fermionic part and by for the bosonic 
one, with the result P] 

aJ^^(two - loops) = -4.4(0.4) ■ 10"^°. (1.3) 

With the same low theoretical and experimental uncertainties, and when com- 
bined with other high precision results, from the LEP experiments and elsewhere, 
would become an excellent probe of physics beyond the Standard Model (extra 
Wr^l gauge bosons, Z' bosons, extra Higgs bosons, SUSY, . . .). For a study of 
the sensitivity of to new physics see for instance 0. 

In the Standard Model, the contributions to fall into three categories: the 
pure electromagnetic (QED) contributions, the electroweak contributions dis- 
cussed above and the hadronic contributions. The QED contributions have been 
calculated and/or estimated up to order (a/vr)^. They give the bulk of the value 
of a^. For an updated value of ci^^° see [|| and references therein. A review of 
the QED calculations is in 0. The result is 

aj™ = 11 658 470.6(0.2) • 10"^° . (1.4) 

The main actual source of theoretical uncertainty still remains in the hadronic 
contributions. The leading hadronic contributions are of two types: the vacuum 
polarization and the light-by-light scattering contributions. The hadronic vacuum 
polarization contributions are the major source of uncertainty at present. Fortu- 
nately, they can be related through a dispersion relation to the experimental ratio 
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R{s) = a{e^e~ hadrons)/(T(e"'"e~ ) |T3- '^^^ planned improvement 

of the experimental determination of R{s) in e.g. BEPC at Beijing, DA$NE at 
Frascati, and VEPP-2M at Novosibirsk will significantly reduce this uncertainty. 
A recent reanalysis of the contribution of the full photon vacuum polarization 
insertion into the electromagnetic vertex of the muon can be found in 
In ^ 



11], using experimental data, they find 



^vac.pol. _ 725.04(15.76) ■ 10' 
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;i.5) 



In |]T2[, using experimental data below 2 GeV^ and accurate QCD calculations 
above this scale, they find 



a 



vac.pol. 



711.34(10.25) • 10 



-10 



1.6) 



An alternative attempt to compute it at lowest order (a/yr)^ within the same 

At order (a/vr)^ 



low energy model we use in this work can be found in | 
there appear other hadronic vacuum polarization contributions that also can be 
expressed as a convolution of R{s) They have been estimated to be [|T5|, [16 



i™" P°'- (higher orders) = -19.9(0.4) ■ 10" 
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Unfortunately, the hadronic light-by-light scattering contribution cannot be 
related to any observable and hence we must rely on a purely theoretical frame- 
work to calculate it. There have been several attempts to calculate this contri- 
bution in the past |ll5| , [T^. There has also been recently some discussion about 
the reliability of this calculation [|17[ |1^. To pin down this contribution is an 



important issue since a quick estimate yields that it could be of the same or- 
der of magnitude as the expected BNL uncertainty. Recently, with the aim of 
reducing as much as possible the theoretical uncertainty from this contribution, 
there have appeared two works, |]T9| and [^, which calculate a^^^^^-^y-^^s^^ , An 
extended version of [19 



IS m 



21 1 . This paper is the detailed version of |20 



The present work is devoted to the calculation of the contributions of the 
hadronic light-by- light scattering to a^. A first simplified version and summary 
of the main results of this paper was presented in [^. A numerical mistake was 
discovered in the first reference of which was corrected in the Erratum. Of 
course, the methods used in 
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and are similar and mainly based on the 
analysis performed in . Nevertheless, we want to discuss the main differences 
in the calculations as well as the reasons why we use the Extended Nambu-Jona- 
Lasinio (ENJL) model for this task. 

The framework we have adopted to calculate the hadronic light-by-light con- 
tribution is an 1/Nc expansion within the ENJL model. In Ref. the 0{Nc) 
leading hadronic contributions were presented. The next-to- leading in the 1/Nc 
expansion effects of the U(l)yi anomaly were also included. There, we took as 

^This number only contains the corrections not included in ( |1.5| ) or ( |1.6| ). 
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a first estimate the result and error for tlie 0{1) in tlie 1/Nc expansion Tr"*" and 
loop contributions from Here we shall come back to all these issues in 
a more detailed fashion. 

The paper is organized as follows. In Section ^ we present the definitions 
related with the hadronic light-by-light scattering contribution to and specify 
the method we have used to calculate it. In Section |^ we explain why we have 
chosen the ENJL model as a good low energy hadronic model. Its main features 
and definitions needed are also presented here. Section ^ concerns the calculation 
of the large Nc contributions to the hadronic light-by- light scattering to a^. This 
section also includes an estimation of the main next-to-leading in 1/Nc effects 
coming from the U(l)^ anomaly. Various checks performed and numerical com- 
parison with other works are also shown here. In Section ^ the next-to-leading 
((9(1)) in 1/Nc contributions coming from charged pion and kaon loops are dis- 
cussed. Here, inspired by the ENJL model 0{Nc) calculation, we will use lowest 
order Chiral Perturbation Theory (CHPT) modulated with vector meson propa- 
gators to calculate them. We discuss issues of gauge and chiral invariance there 
as well. Then in Section ^ we shall gather the numerical results for the contribu- 
tions analyzed in the previous sections. In Section ^we discuss the contributions 
coming from the intermediate (between 0.5 GeV and 4 GeV) and higher energy 
regions to the hadronic light-by-light scattering contribution to a^. In particular, 
how to estimate them and the theoretical error they induce. Finally, Section ^ 
is devoted to a discussion of the results and conclusions. Appendices collecting 
some analytical expressions are also included at the end. 

2 The Method and Definitions 

The amplitude describing the interaction of a momentum p fermion with an 
external electromagnetic field with momentum transferred ps = p — p', can be 
written as 

M = -lelA^uip') 

where m is the fermion mass. The form factor d = —^3(0) is the electric dipole 
moment and fi = -Fi(O) + -^2(0) is the magnetic moment of the fermion in mag- 
netons. In the Born approximation -Fi(O) = 1 and -^2(0) = -^3(0) = 0. In analogy 
with the classical limit, it is usual to define the gyromagnetic ratio g = 2fi and 
the anomalous magnetic moment as a = {g — 2)/2. 

The hadronic light-by-light scattering contributes to a at order (a/vr)^. This 
is a vector four-point function made out of four quark vector currents attached 
to the fermion line with three of its legs coupling to photons in all possible ways 
and the fourth vector leg coupled to the electromagnetic external source. One 
hadronic light-by-light scattering contribution is depicted in Figure The cross- 
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uip), (2.1) 
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Figure 1: Hadronic light-by- light contribution to a^. The bottom line is the 
muon line. The wavy lines are photons and the cross-hatched circle depicts the 
hadronic part. The cross is an external vector source. 

hatched blob is the hadronic part. The momenta in this figure correspond to the 

first permutation of the three vector legs attaching to the fermion. There are five 

more permutations. 

To extract the hadronic light-by-light scattering contribution to the muon 
= {g^ — 2)/2 we shall closely follow Ref. |^ and use the first permutation 

shown in Fig. |l] as an example. The hadronic light-by- light contribution to M. is 

M = leU ^^'^^ ^^'^^ ^ 



(2.2) 



(27r)4 J (27r)4 q^pjplipl - m^) {pi - m^) 
X n^'^"^ {pi,p2,P3) u{p')-fc,{^^ + m)-f^{^r, + m)-fpu{p) 
+ five more permutations , 



where we have explicitly given the contribution from the permutation in Fig. |I[ 
In Eq. (p^) , m is the mass of the muon and n'"'"^(pi,p2,P3) is the four-point 



function 



I E / A / d^ze^(^--+^-^+^-^)(0|r(Kf(0)n'^(a;)V;"(i/)Kf(z)) |0). 



6 , w 

a,b,c,d 



4 



(2.3) 

V/^{x) = Qi [qi{x)'^^qi{x)], qi is a quark of flavour i and Qi its electric charge in 
units of |e|. Summation over colour between brackets is understood. The explicit 
calculation of n'"^"^(pi,p2,P3) in the Standard Model is the subject of the next 
sections. 

The reason for the factor 1/6 in the definition of n^'^"'^(pi,p2?P3) is the fol- 
lowing. This four-point function has six contributions due to Bose symmetry of 
the vector legs, but these permutations are just giving the same six permutations 
we already took into account in Eq. (|2.2|). What we do is to calculate with both 



the six permutations of the four-point function vertices and the six permutations 
of the points where photons connect to the muon line and therefore we have to 
divide by six the four-point function not to make double counting. This will be 
very useful since we can now use both U(l) gauge invariance in the construction 
of n''^°^(pi,P2,P3) (see App. ^ and in the electromagnetic vertices on the muon 
line. 

Because of the U(l) gauge invariance, the sum over all six permutations con- 
tributing to I\.P'^°'^{pi,p2iP'i) has to be UV finite in renormalizable theories as 
the Standard Model. This is true for the sum of all the permutations of a given 
class of contributions to n''''"^(pi,p2,P3)- However, each single permutation can 
be divergent. This is the case for the one fermion loop contribution where some 
permutations are UV logarithmically divergent. In order not to rely on numerical 
dangerous cancellations, we write n^''^°'(pi,p2,P3) in a form where each permu- 
tation is explicitly UV convergent. Following Ref. we use the U(l) gauge 
covariance condition p3^n''^"^(pi,p25P3) = to obtain 

n--(P.P.P3)^-P3.^5q!««^. (2.4) 



The presence of the extra derivative with respect to makes the rhs of (|2.4| ) 
explicitly UV finite. Then M. can be rewritten as 

M = -\e\AxP3(3u{p')M^''{ps)u{p), (2.5) 

where, for the permutation in Eq. ( |2.2| ), 



(27r)4 J (27r)4 q^ pj pl{pl - m^) [pj - m^) 
5nP-/3(p,,p2,p3)" 



6p: 



'3A 



It can be shown, see [22|, that the contribution from the light-by-light scattering 



to can be written as 



^light-by-iight ^ tr[(^+m)M^'^(0)(^ + m)[7A,7/3]]. (2.7) 

4oT7l 
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The eight dimensional integral in Eq. ( p.6| ) of the two loops on muon momenta 
can be reduced to a five dimensional one, two moduli and three angles, using the 
symmetries of the system. The integration over these variables has been done in 
Euclidean space. The integral in Eq. (|2.3| ) brings in other integration parameters, 
the latter integral we also perform in Euclidean space. 

The momenta flowing through the three photon legs attached to the muon 
line run from zero up to infinity, covering both the perturbative and non per- 
turbative regimes of QCD. These two different regimes are naturally separated 
by the scale of the spontaneous symmetry breaking ~ 1 GeV. In the region 
between A-^ and say (4~5) GeV, the strong interaction contributions have to 
match the perturbative QCD predictions in terms of quarks and gluons. Hence, 
rigorously, one should calculate the low energy contribution to n''^°^(pi,p2 5P3) 
and match this result with a perturbative QCD calculation of the high energy 
contribution. Here this is technically rather involved because of configurations 
with both high and low energy photon momenta. If the scale determining the 
bulk of the contributions to a''g'^*~'^y~''g'^* were around the muon mass we could 
attempt to make a pure low energy calculation that would saturate it. Were the 
contributions from the high energy perturbative region not negligible we would 
need a more sophisticated model suitable both in the low and intermediate energy 
regions. We have investigated this issue by putting relevant Euclidean ultraviolet 
cut-offs, labelled fi from now on, on the moduli of the momenta attaching to the 
muon line in Eq. ( p.6|) . 

Since, as said before, the momenta in the photon legs can run up to infinity any 
accurate calculation should incorporate the full external momenta dependence of 
llP'"^l^(^p^^P2,p^). This is clearly beyond the reach of the present state of the art 
of CHPT. Alternatively, one can rely on a good low energy hadronic model. We 
have chosen the Extended Nambu-Jona-Lasinio model. Reasons for that choice, 
definitions and main features of the model are in the next section. By using this 
model, we calculate the low energy contributions to a|^gi^*-t'y-iisht^ then study 
the saturation of a^^^^^-^^-^^st^^ by the physics at scales below or around A, where 
A in our case is the physical cut-off of the ENJL model. Although it will turn 
out that the contributions from intermediate and high energy regions are not 
negligible, we shall be able to give a conservative estimate for them. This will be 
explained in Section |^. 



3 The ENJL Model 

For recent comprehensive reviews on the NJL [p3| and the ENJL models [E| 



see Refs. Here, we will only summarize the main features, notation and 

reasons why we have chosen this model. More details and some motivations on 
the version of ENJL we are using can be found in ^7\, ^ . 

The good features of the ENJL model and its suitability for an accurate 
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calculation of aj^g'^*^^^^''^'^* were realized in Ref. [|13]. In fact, in this reference the 
hadronic vacuum polarization to lowest order (a/vr)^ was calculated within the 
same ENJL model we use here, obtaining a good agreement (~ 15%) with the 
phenomenological result in Ref. |T^. This accuracy, enough for our purposes, 
is the maximum we can expect from our calculation. One of the conclusions 
of reference was that the hadronic vacuum polarization to saturates for 
energies around 1.5 GeV, which is still a reasonable scale where the model could 
be applied without introducing too much uncertainty. We shall see that in the 
light-by-light case the contributions from higher energies are not negligible. More 
comments on this are in Section ^ 

3.1 The Model and Determination of Parameters 

The QCD Lagrangian is given by 

-^QCD = 1 {^1^ (d,, - iv^, - ia^-f5 - iG^) - {M + s - ijry^)} q ■ (3.1) 

Here summation over colour degrees of freedom is understood and we have used 
the following short-hand notation: q = (u, d, ; is the gluon field in the 
fundamental SU(Nc) (Nc=number of colours) representation; G^u is the gluon 
field strength tensor in the adjoint SU(Nc) representation; v^, a^, s and p are 
external vector, axial-vector, scalar and pseudoscalar field matrix sources; Ai 
is the quark-mass matrix. The ENJL model we are using corresponds to the 
following Lagrangian 



Cenjl = C^ciy + '^9sT.{^yL){^L<tR) 



-9v E [{tL^qi) {qii.ql) + (^ - R)] ■ (3-2) 



Here i,j are flavour indices, \E'k,l = (1/2) (1 ± 75) ^ and 



ill 



^Gv{A) _ 47r2G's(A) 



The couplings Gs{A.) and Gv{A.) are dimensionless and 0{1) in the l/N^ expan- 
sion and summation over colours between brackets in (|3.2| ) is understood. 

The Lagrangian Cqqj) includes only low frequency (less than A) modes of 
quark and gluon fields. The low frequency modes of the gluon fields can be 
assumed to be fully absorbed in the coefficients of the local operators or alterna- 
tively described by vacuum expectation values of gluonic operators. So at this 
level we have two different pictures of this model. One is where we have in- 
tegrated out all the gluon degrees of freedom and then expanded the resulting 
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effective action in a set of local operators with quark fields keeping only the first 
non-trivial terms in the expansion. The other picture is obtained when we only 
integrate out the short distance part of gluons and quarks. We then again expand 
the resulting effective action in terms of low energy local operators with gluon 
and quark fields. This is described in |^ ^ and the best fits there correspond 



to the first alternative. Therefore, in the present work we will use (|3.2|) with all 
gluonic degrees of freedom integrated out. 

This model has three parameters plus the current light quark masses. The 
first three parameters are Gs, Gy and the physical cut-off A of the regularization 
that we chose to be proper-time. Although this regulator breaks in general the 
Ward identities we impose them by adding the necessary counterterms (both in 



the anomalous and in the non anomalous sectors). The light quark masses 
in Ai are fixed in order to obtain the physical pion and kaon masses in the poles 
of the pseudoscalar two-point functions |]2B|. The values of the other parameters 
are fixed from the results of the fit to low energy effective chiral Lagrangians 
obtained in They are Gs ^ 1.216, Gy ^ 1-263, and A = 1.16 GeV from 



Fit 1 in that reference. Solving the gap equation, we then obtain the constituent 
quark masses: M„ = = 275 MeV and = 427 MeV. 

3.2 Advantages and Disadvantages of the Model 

The ENJL model is a very economical model that captures in a simple fashion 
a lot of the observed low and intermediate energy phenomenology. It has also a 
few theoretical advantages. 



1. The model in Eq. ( p.2| ) has the same symmetry structure as the QCD 
action at leading order in l/Nc [0. Notice that the U(1)a problem is 
absent at this order [^]. (For explicit symmetry properties under SU(3)l 



X SU(3)j? of the fields in this model see reference l^^-) chiral limit 

and for Gs > 1 this model breaks chiral symmetry spontaneously via the 
expectation value of the scalar quark- ant iquark one-point function (quark 
condensate) . 

2. It has very few free parameters. These are unambiguously determined 
from low energy physics involving only pseudo-Goldstone bosons degrees 



of freedom[27 . 



It only contains constituent quarks. Therefore, all the contributions to a 
given process (in particular to a^) are uniquely distinguished using only 
constituent quark diagrams. Within this model there is thus no possible 
double counting. In particular, the constituent quark-loop contribution 
and what would be the equivalent of the meson loop contributions in this 
model, are of different order in the 1/A^c counting |[T^. This shows clearly 
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(a) (b) 

Figure 2: The graphs contributing to the two point-functions in the large Nc 
hmit. a) The class of all strings of constituent quark loops. The four-fermion 
vertices are those in Eq. ( p.2|) . The crosses at both ends are the insertion of the 
external sources, b) The one-loop case. 



that there is no double counting here. As described in ^Tj this model 
includes the constituent-quark loop model as a specific limit. 

4. Two-point functions are given by the general graph depicted in Fig. |^. The 
resummation of strings of constituent quark bubbles for two-point functions, 
automatically produces a pole in all main spin-isospin channels within the 
purely constituent quark picture. This is qualitatively the same as the 
observed hadronic spectrum. 

5. It provides a reasonable description of vector and axial vector meson phe- 
nomenology. 

6. Some of the short distance behaviour is even the same as in QCD. For 
instance, the Weinberg sum rules [BBl are satisfied. This allows in some 



cases for a good matching with the short distance behaviour, see below for 
an example. 

7. The major drawback of the ENJL model is the lack of a confinement mech- 
anism. Although one can always introduce an ad-hoc confining potential 
doing the job. We smear the consequences of this drawback by working 
with internal and external momenta always Euclidean. 

The techniques used here together with more phenomenological issues are 



treated in Refs. BE, 32, 33] and reviewed in pq]. Some applications to other 



non-leptonic matrix elements can be found in 0] and The general conclu- 



sion is that within its limitations the ENJL-type models do capture a reasonable 
amount of the expected physics from QCD, its symmetries, their spontaneous 
breakdown and even some of its short distance information. 

As mentioned before, the Weinberg sum rules are satisfied. This is a 
very important point and is another of the reasons why we have chosen this 
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model. These relations are needed to obtain good matching between the low- 
energy behaviour and the high-energy one. As an example, they are essential 
for the convergence of the hadronic contribution to the electromagnetic tt"*" — 7r° 
mass difference Models to introduce vector fields like the Hidden Gauge 

Symmetry (HGS) |^ do not have this good intermediate behaviour for some 
choices of the parameters. The choice of parameters in the HGS model used in 
1^ to calculate is affected by this problem. For instance, the contribution to 



the above mass difference in the HGS model up to a cut-off z/, is [39 



3 a 

4 TT 



'1 



a)i/^ -aMy log [My/ 1 



(3.4) 



where a is a parameter of the HGS model. This obviously diverges badly for 
a = 2 which is the value chosen in [|19|. The same problem is avoided in the 
ENJL model ISl. 



3.3 The Present Usage. 

We will use the ENJL model as a model to fairly describe in the large Nc limit 
strong interactions between the lowest-lying mesons and, if needed, external 
sources. This is a tree-level loop model with an explicit cut-off regularization 
for one loop parts. What we mean by a tree-level loop model is the following: 
a general set of external sources is connected via full chains, like the one de- 
picted in Fig. 0, to one-loop diagrams which are also glued through full chains 
or four-fermion ENJL vertices. These are the leading contributions in l/Nc- It is 
at this level that the hadronic properties of this model have been tested. To go 
beyond this level one would have to include other operators not suppressed at the 
next-to- leading order in 1/Nc in the ENJL Lagrangian. At that level one also en- 
counters the problem of regularizing overlapping divergences in the model. This 
is the reason why we shall only apply the ENJL model to calculate the low-energy 
large contributions to ajs^^^-^y-^^s'^*. 

Since one of the issues that motivated this calculation was the apparent not 
fulfilling of Ward identities in previous calculations (see comment in |1^), we want 



to emphasize here that this model possesses chiral symmetry and the necessary 
counterterms are added so that n— point Green functions fulfil both anomalous 



and non anomalous Ward identities 29 . For instance, the calculation in 19 



assumes ordinary VMD for the anomalous sector. It was shown in that this 
VMD breaks the anomalous Ward identities and one needs for subtractions to 
restore them. This is particularly important for the flavour anomaly contribution 
to the hadronic light-by-light scattering. We want to point out also that in 
the ENJL model we are using, both anomalous and non anomalous sectors are 
described by the same set of parameters. This is not the case for HGS models 
where consistency between parameters in both sectors is not obvious. 
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4 Low Energy Large Nc Contributions 

In this section we discuss the low energy contributions that appear at large Nc 
within the framework of the ENJL model |]T3. For a general contribution, as 



can be seen in Eqs. ( p.6|) and ( |2.7| ), we have to compute the derivative of the 
generalized four-point function n^'^"^(pi,p2,P3) in Eq- (|2.3| ) with respect to p^x 
at ps = 0. The Lorentz structure of this four-point function and some other 
technical aspects of its calculation are in Appendix ^. Since we are dealing 
with the low energy contributions to we shall only consider the lightest quark 
flavours: up, down and strange. Contributions from heavier flavours are discussed 
in Section |^ In the ENJL model there are two classes of contributions to the 
generalized four-point function n^^"^(pi,p2!P3)- The first one is a pure four- 
point function (see Fig. 0(a)). The second one, which we call three-point-like 
function contributions, can be regarded as two three-point functions glued with 
one full propagator (see Fig. ||(b)). Within this framework, it clearly appears 
that these two classes of contributions to the generalized four-point function 
YlP^'^l^(^p^^P2^p^) are different and therefore should be summed up [|l^. The pure 



four-point function corresponds to the so-called quark loop contribution and the 
three-point-like function contributions to the meson pole exchange in the language 
of Ref. [|T^. We remind that in our framework both classes of contributions to 
-QP!/a/3^p^ g^j^g calculated to all orders in the CHPT expansion. 

4.1 Pure Four-Point Function Contribution 

This contribution is diagrammatically represented in Fig. ^(a). In the CHPT 
expansion the lowest order contribution is order p®, thus potentially sensitive to 
the high energy region. The momenta assignment shown in the figure is the one 
corresponding to the first permutation. Whenever we give explicit expressions 
these are the ones corresponding to this permutation. The other five permutations 
are done analogously. 

The pure four-point function is formed by a one-loop constituent quark vector 
four-point function H'^""'^ {pi,p2,P3) with the same Lorentz and flavour structure 



as the four-point function in Eq. (|2.3|) with the three vector legs attaching to 
the muon line dressed with full vector propagators. These are the cross-hatched 
blobs in Fig. |](b). The circled crosses vertices in the figure are where the photon 
lines connect. In the large Nc limit there is either a constituent up, down or 
strange quark running in the loop. For a given quark flavour, this contribution 
to n'"^°^(pi,p2 5P3) can be written using the methods of pSf as, 

n''-^(pi,P2,P3) = TiaUpl,P2,P3) V^'^^'^^ipi^p^^p,) + ■■■ (4.1) 

with 

V M^{-q^)-q^ / \ M^{-pi)-pi I 
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Figure 3: The two classes of hadronic light-by-light contributions to at leading 
0{Nc). (a) The four-point functions class, (b) The product of two three-point 
functions class. The dots are ENJL vertices. The circled crossed vertices are 
where photons connect. The cross-hatched loops are full two-point functions and 
the lines are constituent quark propagators. 
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X 



(4.2) 

Dots stand for other possible contributions to 'n.P'^°'^{pi^p2.,Pz)- The indices 
a, 6, c, d, e, / stand in the remainder for Lorentz indices. The function My{—p^) 
corresponds to either mm, dd or ss flavour structure and can be found in Ref. [Q. 
Notice that the full vector propagators in the ENJL model have the same form as 
in VMD models but with a momentum dependent "vector mass" . See discussion 
in about meson dominance in ENJL models. 

After summing over all the six possible permutations of the external momenta 
only the terms proportional to g^j^y tensor survive because of the U(l) gauge 
invariance. This leads then to the phenomenological VMD rule of replacing the 
propagator of a photon with momentum p by 



p2 M^-p2 



(4.3) 



Thus implying that this rule preserves the Ward identities after summing over 
all the permutations. 

We compute the one-constituent quark loop following the analysis in Ap- 
pendix ^ and using a proper time regulator with a physical cut-off A. 

The contribution of 11 (pi,P2,P3) in eq. ( ^?T| ) to can be decomposed in 
terms of 32 independent amplitudes (see App. 0). After integrating over the 
momentum running in the loop, these amplitudes will be integrals over the three 
Feynman parameters introduced by the standard procedure of reducing the four 
internal propagators to a power of one propagator. Moreover, after taking the 
derivative with respect to p^x at p3 = one of the Feynman parameters can be 
integrated out analytically. The general form of each one of the 32 functions in 
n (^1,^25^3) for a constituent quark of flavour i contributing to a|^s'^t-t'y~iis'^t 
is thus 

where 

M2(x, y) = M^ - x(l - x)(pi + p2f - y{l - y)pl + 2xy{pi + P2) ■ P2, (4.5) 
and 

^2(e) = (l + e)e-^ (4.6) 

Here Mj is the constituent quark mass with flavour i. The Feynman parameters x 
and y, together with the other five degrees of freedom for the two external muon 
loops in Eq. ( p.6| ), produce a seven dimensional integral that we perform using 
the Monte Carlo routine VEGAS. The numerical results of this contribution will 
be discussed in Section ^ 



13 



Cut-off 


X 10'^ 


X lO'^ 


ttf, X 10^ 


jj, 


Electron 


Muon 


Constituent Quark 


(GeV) 


Loop 


Loop 


Loop 


0.5 


2.41(8) 


2.41(3) 


0.395(4) 


0.7 


2.60(10) 


3.09(7) 


0.705(9) 


1.0 


2.59(7) 


3.76(9) 


1.10(2) 


2.0 


2.60(6) 


4.54(9) 


1.81(5) 


4.0 


2.75(9) 


4.60(11) 


2.27(7) 


8.0 


2.57(6) 


4.84(13) 


2.58(7) 


Known Results 


2.6252(4) 


4.65 


2.37(16) 



Table 1: The results from a pure fermion loop compared with known numerical 
and analytical results. The known result for the electron loop has been partially 
analytically evaluated in pO|, while the quoted values for the muon loop and the 



constituent quark loop are respectively analytical and numerical results of |16] 



As a check we have reproduced the results for the constituent quark and muon 
loops in |T^ and the electron loop result in [^. This is the calculation using 
([4.1|) with V"^cdpi/a/3 _ gap gbu gca gdp ^ have also re-scaled the relevant masses. 
For the constituent quark mass we have used Mq = 300 MeV. The results are in 
Table The numbers in brackets are the VEGAS quoted errors. The last row 
gives the known results. Within the numerical uncertainty of our calculation the 
agreement is good. It should be remarked that there are cancellations present at 
high energy. These contributions make the estimate numerically more uncertain 
for higher values of the cut-off /x. As can be seen the electron loop has essentially 
reached its final value at a cut-off of 500 MeV. The same is not true for the muon 
loop or the constituent quark loop. The muon loop reaches its asymptotic value 
essentially at a cut-off of about 2 GeV. The constituent quark loop reaches its 
asymptotic value above 4 GeV. We have also calculated the same quantities for 
a few higher cut-offs and observed their stabilization within the numerical error. 
These sample calculations show that a proper study of the cut-off behaviour 
for the more difficult to estimate complete hadronic contribution is definitely 
necessary. 



4.2 Three-Point-Like Function Contributions 

This class of contributions is diagrammatically represented in Fig. ^(b). There 
are two permutations of the vector legs for each of the two three-point functions. 
In addition there are three possible sets of two momenta out of the four external 
vector legs momenta. This makes twelve possible permutations to be considered 
for three-point-like function contributions. The momenta shown in the figure 
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are the ones corresponding to the first permutation. Whenever we give exphcit 
expressions these are the ones corresponding to this permutation. The other 
eleven permutations are done analogously. 

In this case we have two one-loop three-point functions with two vector legs 
each one glued with a full two-point function that can be either pseudoscalar, 
scalar, mixed pseudoscalar-axial-vector, or axial-vector. For intermediate vec- 
tor two-point functions the result is zero because of Furry's theorem. Three of 
the vector legs here are then attached to the muon line with dressed full vector 
propagators just as in the case of the pure four-point function discussed in the pre- 
vious section. Again, in the large Nc limit, there is either an up, down or strange 
constituent quark running in the loop. Technically we have used two different ap- 
proaches to calculate this type of contributions. One is using the Ward identities 
for the four-point function Il'""^^{pi,p2,P3). In this way one has to determine the 
32 amplitudes of n''^"^(pi,p2,P3) contributing to ajj^'^*"^^"''^'^*, see Appendix ^ 
The other way is constructing explicitly the complete Il'"^°'^{pi,p2,P3) starting 
from the three- and two-point functions. Here one relies on the Ward identities 
for three-point functions. As a check we verified that both ways agree exactly. 

In what follows we study each type of exchange (scalar, pseudoscalar and 
axial-vector) separately. 

As an additional numerical check we have also calculated the pion exchange 
contribution with vector meson dominance of the rho meson in the vector legs 
and when the anomalous vertices are from the order p'^ Wess-Zumino effective 
action. Our result agrees exactly with that in Eq. (4.1) of |^ which quotes 
-55.60(3) ■ 10-^\ 

The two-point function involves the integration over one Feynman parame- 
ter and each one of the two three-point functions the integration over two more 
Feynman parameters. These integrals have been evaluated using Gaussian inte- 
gration. To obtain a^, one has to convolute these three-point-like contributions 
with the five dimensional space integral of the external two muon-loops in Eq. 
( |2.6| ) which has been performed using the Monte Carlo routine VEGAS. The 
numerical results for these contributions will be discussed in Section ||. 

4.2.1 Scalar Exchange 

For the scalar exchange, the lowest order contribution in the CHPT expansion is 
order p^. This contribution to n''^"^(pi,p2,P3) can be written as 



W^-P{p„P2,Ps) = U:r{Pur)gs{l + gsIl'{r))u:riP2,P3)V-'^'''''^Pi,P2,P3) 

+ ■■■, 

(4.7) 

where V^'^^'^p'^'^I^ (pi , P2 , Ps ) has been defined in Eq. ( ^4.2| ) and r = p2 + ps- The 
full two-point function Il^{p) can be found in Ref. and the one- loop three- 
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point functions 11^^ {p, q) and 11^^ (p, q) are in Appendix 0. Eq. ( [4.7|) can be 
understood as follows: gsi)- + gs^^ij)) is the scalar propagator connecting the 
two vertices 11^^ (^2,^3) and 11^^ {pi,r). The whole diagram is then glued to 
the external photon lines with the vector propagators given by V'^^cdpj/a/?^ 



4.2.2 Pseudoscalar Exchange 

For the pseudoscalar exchange the lowest order contribution in the CHPT ex- 
pansion is order p^. This, together with the fact that it involves two flavour 
anomalous vertices points out that this contribution could be the leading one. 
Considering as part of the pseudoscalar exchange all those terms proportional 
to a pseudoscalar propagator, this contribution also includes the pseudoscalar- 
axial-vector mixed terms. Its expression is given by 



n^^"^(pi,P2,P3) = [llir{Pi.r){l + gsTl''{r))TlT{p2,p,) 



X 



-9v^lb^{Pi. r)lf'"{r)U:l^^{p2, ps) 
-9v n^fe^^(pi, r)n^^(-r)nf7^(p2, P3) 

^75V"'^'^""^(Pl,P2,P3) + --- 



(4.8) 



where V"''^'^''''"^ (Pi , P2 , P3 ) has been defined in Eq. 



and r = P2 + ps- The 



full two-point functions H (p) and n^(p) and the one- loop three-point functions 

51) in Appendix 



nfj^(p, q) and U^^^{p, q) can be found in Ref. [Q. See Eq. 



|B|for the explicit expression. The mixed two-point function n^(j9) can be written 
as PI, 

Ill{p)^ip,Ill,{p'). (4.9) 

Since both the VVA and the AVV one-loop three-point functions are multiplied 
by r^, we use the one- loop anomalous Ward identity in Eq. (4.24) of pSf and the 
prescription given in to rewrite Eq. ( |4.8| ) as 

n''^'^^(pi,P2,P3) 



gslla, {Pu r) (1 + ^?5n^(r) - ^vM.UU 



'-ab 

ttPVV 
'-cd 



xn:;'(p2,P3)v''^^'^'"^"^(pi,P2,P3) 



+ 2MQgsgvKA 
fg'^mU-q'] 



n 



VVP 
ab 



q"-qp 



{pi,r) {n, 

fg^^M, 



t'vi 



PVV 

cd 

-pI) 



(P2,P3)' 

-pIpi 



9 9 



MU-l'')-<l'' 



,Mi=MQ 



mU-pD-pI 



+ {^vvpiPur) 



^PVviP2,P3) 



gapgtu 



gcagdp 



' g'^Ml{-pl)-p'^^p^ \ ( g''^M^{-pl)-pipl \ 



m^-pD-pI j 



mU-pD-pI J 



+ 



(4.10) 
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where Mj is the constituent quark mass for the quark with flavour i and Mq 
its value in the chiral limit. The first two lines above come from applying naive 
Ward identities to the axial-vector leg, while the remaining contributions are 
the subtractions needed to fulfil the anomalous Ward identities. In addition, 
(p, q) contains subtractions also determined by the anomalous Ward iden- 
tities p8|. See Eq. ( p.6| ) in Appendix 0. The expression (|4.10| ) also shows that 
the pseudoscalar exchange contribution always contains at least one vector meson 
propagator. Due to this fact and the presence of the subtractions in ( [4.10|) , the 
77*^7*7* vertex goes to a constant when the vector legs' Euclidean momenta are 
very large. Therefore, although this contribution when summed over all possible 
permutations is convergent by itself because of gauge invariance (see Section 
the subtraction terms make it very slowly convergent. We know that in QCD the 
7r''7*7* vertex goes like 1/Q^ at large Euclidean momentum [^. This behaviour 
is also supported by the measured 77*^77* form factor in the Euclidean region at 
CELLO |42| and CLEO-II |43| detectors. This indicates again that although this 
model gives the right contribution for energies below or around A, it breaks down 
above. We shall estimate the intermediate and high energy region contributions 
for the pseudoscalar exchange in Section ^. 

Although the present section is devoted to the large Nc contributions, it is 
worth to discuss the main 1/Nc corrections to the pseudoscalar exchange at this 
point. These are the effects of the U(l)^ anomaly and were already included 
in the Erratum in Ref. In the chiral limit and in the large limit there 

are nine pseudo-Goldstone bosons [^: tt^, 7r°, K^, K^, , rjg and r/i. Under 
flavour SU(3) they transform as a nonet multiplet. However nonet symmetry is 
broken by l/Nc effects due to the U(1)a anomaly. These effects cause the isospin 
zero mass eigenstates to become the rj and 77' states. They also increase the mass 
of the 7]' meson to 958 MeV. These l/N^ corrections are thus quite relevant to the 
pseudoscalar exchange. We have taken them into account by using the physical 
7r°, 7] and f]' mass eigenstates as propagating states. This already gives the bulk 
of the effects of the U(1)a anomaly. Higher order corrections are negligible and 
within the quoted error. The results of using either uu, dd, and ss basis of states 
for the large Nc limit or the physical 7r°, rj and 1]' basis are given in Section 0. 



4.2.3 Axial- Vector Exchange 

For the axial- vector exchange, the lowest order contribution in the CHPT expan- 
sion is order p^. This contribution to Il'""^^{pi,p2,P3) can be written as 

n^'^"^(pi,P2,P3) = -9vTl2^iPi,r){ff-gvIl':iir))TlZ''iP2,P3) 

y<V^''''"'"''iPi,P2,P3) + --- (4.11) 

where V^bcdpuafS ^p^^p^) 

has been defined in Eq. (|4.2| ) and r = P2+P3- The full 
axial- vector two-point function Il'^{p) can be found in Ref. and the one-loop 
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three-point functions Il^^J^{p,q) and ]l^JJ(p,q) are in Appendix 0. 

Using the one-loop anomalous Ward identity in Eq. (4.24) of and the 



prescription given in for the terms in ( 4.11|) where multiplies either the 



VVA or the AVV three-point functions, we can rewrite ( [4. Ill ) as 

r / ^"^M^(-g^)-gV \ ( g''MU-pl)-plp'i \ 
MjiTiPu r) - MqH- r) ^^2^^, } 

r / (?^°M^(-p^)-pip^ \ / ^'^^M^(-p^)-pfpg \ 

U M2(-pi)-pi ; 1, M2(-pi)-pi ; 



X 



+ (4.12) 



The expression for the two-point axial amplitudes 11^^ (—r^) and 11^^ (—r^) can 
be found in |2^. The kinematical pole of the transverse part at = disap- 
pears in the combinations 1 + gyr'^Ii^^ {—r'^) and 11^^ (—r^) -|-n^^(— r^). The last 
combination contains a pseudoscalar propagator which is often included in the 
pseudoscalar exchange contribution. 



5 Low Energy Next-to-Leading in 1/Nc Contri- 
butions 

In this section we discuss the low energy next-to-leading in l/N^ contributions 
to aj^^'^*~''^~^'^^*. In the previous sections we have discussed the low energy large 
Nc contributions in the context of the ENJL model. We included the effects of 
the U(l)^ anomaly which are next-to-leading in l/Nf. as well. Here we address 
the calculation of the other 0{1) corrections in the l/N^ expansion. In the 
language of the ENJL model, these corrections are given by loops of strings of 
quark bubbles. They contain one closed loop of a string of bubbles like the one 
in Fig. @(a). This loop is then connected in all possible ways to the photons 
using strings of bubbles. In a mesonic picture they correspond to one meson loop 
contributions. This meson loop can then substitute any one-loop constituent 
quark (bubble) in diagrams in Figs, ^(a) and 0(b). We end up with two classes 
of 0(1) contributions. The 1/Nc next-to-leading corrections to the full two- 
point functions, i.e. to meson masses and couplings and the l/N^. next-to-leading 
corrections to the three- and four-point one-loop functions, i.e. to vertices. In 
fact, the major l/N^ correction to pseudoscalar two-point functions, namely the 
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U(l)yi anomaly, was already estimated in Section ^. The other 1/Nc corrections 
to two-point functions should be small since the phenomenological analysis in 
Refs. |3^ and ||28[ fits very well. We thus expect these to be already included 
in the error of the model for the large Nc results. Therefore, we shall only consider 
here the 1/Nc corrections to the vertices. 

Unfortunately, at present, these type of contributions cannot be fully treated 
in the ENJL model. Some of the reasons were given in Section ^. In its present 
form the ENJL model we are using is just well defined in the large N^. limit. How- 
ever, the fact that these type of contributions are of different large counting 
with respect to the ones treated in Section ||, permits the separate treatment we 
follow below for them. 

Four-point functions can be also calculated at very low energy within CHPT 



1^. In this regime the relevant degrees of freedom are the lowest pseudoscalar 
mesons, while vector, axial- vector and scalar resonances have been integrated out. 
Their effects are included in the couplings of the CHPT Lagrangian CHPT 
becomes then a good tool to study strong interactions of the lowest pseudoscalar 
mesons with external sources. In this framework we need both 7*P+P^ and 
ry*>^*p+p- vertices, where P is pion or kaon.| The first vertex is well known 
phenomenologically and VMD models give a very good description of it. On the 
contrary, not much is known phenomenologically about the second one. This fact 
induces a large model dependence since one can construct many models satisfying 
the relevant Ward identities and with different degrees of VMD. One can use for 



instance the HQS model as in [11^ where there is no complete VMD for the 
^*^*p+p- vertex. We shall discuss this contribution in a complete VMD model 
both for '~f*P~^P~ and 7*7*P"'"P~ vertices. This is done inspired by the form of 
the 0{Nc) contributions in the ENJL model. 

In previous sections we have seen that the most general expression for the 
four-point function H''^"^(pi,p2,P3) at 0{Nc) in the ENJL model contains the 
one-loop four- or three-point-like functions, which give the lowest order in the 
CHPT counting, multiplied by the ENJL vector meson propagators in Eq. ( [4.2| ). 
Inspired by this behaviour we saturate the 0{1) contribution by one loop of 
charged pion or kaon mesons using lowest order CHPT photon-pion vertices, i.e. 
(9(j9^), multiplied by the ENJL vector meson propagators in ( |4.2| ). The main 
difference with a full ENJL calculation here is that we substitute the momentum 
dependent pion mass and coupling by their experimental values. 

To understand the sensitivity to the momenta dependence of the vector me- 
son propagators we have numerically studied the difference between the case with 
vector propagators containing a constant vector mass and the case with a mo- 
mentum dependent one. Each choice gives different high energy behaviour of the 
photon-pion vertices. 



■^Contributions from vertices with more photons start at higher order in CHPT and are 
therefore suppressed with respect to the contributions of 7*P+P^ and 7*7*P+P^. 
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For the pseudo-Goldstone bosons loops contribution to n''^°^(pi,p2,P3), the 
lowest order in the CHPT counting is order p^. This also points out that this 
contribution could be dominated by lower energy regions than the other contri- 
butions analyzed in previous sections which started at order p^. Since we are 
including as propagating states only the lowest pseudoscalar mesons and the rho 
vector, our approach will be only valid for energies below or around 1 GeV. Above 
this energy axial-vector mesons become also dynamical states. The saturation of 
this contribution from physics at scales around 1 GeV can only be confirmed a 
posteriori. The numerical results concerning this contribution are presented in 
Section ||. 

We now proceed to analyze in more detail the two types of possible contribu- 
tions: the pure four-point function and three-point-like function contributions. 
In particular, the three-point-like contributions can be seen as the diagram in 
Figure ^(b) where one or both one-constituent-quark loop three-point functions 
are replaced with charged pion or kaon loops. Due to parity, the intermediate 
two-point function glueing the two three-point functions can be either vector or 
scalar. The vector contribution is again zero because of Furry's theorem (this can 
be better verified in the ENJL inspired diagrams where the vector legs couple to 
fermion lines). The scalar contribution we expect to be very much suppressed like 
in the 0{Nc) case (see numerical results in Section |]). The important point here 
is that there are no anomalous contributions since now we have mesons running 
in the three-point functions. This makes this contribution to be in the range of 
the expected CHPT counting and not anomalously large as we obtained for the 
C(iVc) pseudoscalar exchange. 

Therefore we only estimate the dominant pure four-point function contribu- 
tion to n'"^"^(pi,p2;P3) in Eq- ( p.3| ). This can be seen as the diagram in Fig. 
11(a) where the one-loop constituent quark four-point function is now a loop of 
pseudoscalar mesons. It can be written as 



where ^abcdij>ii'P2,'Pz) is the 0{1) in IjN^ contribution from charged pion and 
kaon loops using lowest order in CHPT photon-pion vertices. For that, we com- 
pute at lowest order in CHPT, the quark vector current appearing in the definition 



where the subscript aa means that we take the aa component in flavour space. 
Due to the covariant derivative D^P = d^P — i\e\[A^,P], the term above gives 
rise both to P~^P~'j* and P^P~'~f*'~f* vertices. Of course, the full contribution to 
]Jp^°'i^(^p^^P2^p^) at a given order in CHPT, in this case 0{p'^), has to be gauge co- 
variant. We construct the one-loop four-point function Ilabcd{Pi,P2,P3) following 
again the analysis in Appendix ^, where we need to determine the 32 independent 



W'"'^{PUP2,P3) = UaUPuP2,P3) V'^''^''''-^ (p^ P2, Ps) + " " " 



(5.1) 



of np-/5(pi,p2,P3) H 




(5.2) 
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amplitudes that contribute to the ajjg'^*"^^"''^'^*. These amphtudes have always 
the Lorentz indices saturated by external momenta indices. Notice that vertices 
with one photon are proportional to external momenta, while vertices with two 
photons are proportional to g^u tensors. Thus we only need to compute the UV 
convergent amplitudes and reconstruct the contribution from the two-photons- 
two-mesons vertices by using gauge invariance. The full meson one-loop with 
order p"^ vertices explicitly satisfies the chiral and U(l) gauge Ward identities. In 
fact the amplitudes one gets for the 32 independent functions are very similar to 
the ones found in the one-loop constituent quark amplitudes in Eq. (^4.4| ). The 
general form we obtain for them is 

(Qu-Q,)' ['dxT^dy (5.3) 



167r2 Jo Jo M\x,y) 

with 

M'^{x,y) = mp- x{l - x){pi+p2f - 2/(1 - y)pl + 2xy{pi+p2) ■ P2 (5.4) 
and mp is the mass of the pseudoscalar meson P. The one- loop meson pure four- 



point function is multiplied with the propagator in ( [4.2| ) to give the full four-point 
function of Eq. ( ^.1|) . 

If the one-loop four-point function satisfies the Ward identities, as it does, 
also the full four-point function ( ^.1|) will satisfy them since 



q, (g^^Ml - q^q^) = q'' {m^ - q') , (5.5) 

independently of whether My is momentum dependent or not. As we already 
said for the one-loop constituent quark loop contribution, after summing over all 
possible permutations of the three vector legs, only the terms proportional to gfj,u 
in ( [4.2| ) survive since the meson one-loop four-point function satisfies the Ward 
identities. This leaves the phenomenologically VMD rule in Eq. ( [4.3|) to work 
here too^. This eliminates the worries about the fulfilling of chiral symmetry 



when using this substitution pR IS 



That this procedure is fully chiral and U(l) gauge invariant can be seen simply 
by constructing a Lagrangian with full electromagnetic gauge and chiral invari- 
ance and that has complete VMD for both 7P+P~ and '~f'~fP^P~ vertices. The 
following Lagrangian contains couplings of pions and photons to all orders in 
external momenta and reproduces the full VMD amplitude without inducing any 
extra vertices of photons and pseudoscalar mesons. 



%MD = ^tr {v,UV^U^ + xU^ + Ux 



Hn fact, the complete VMD we are using is identical to the so-called naive VMD model 
in Ref. |l^ which therefore does not break either the chiral Ward identities or the U(l) 
electromagnetic covariance. 
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n>0 

+ E ^nn^ti [{V^"V,L^^'') [/t (V^^V-R^,) U] . (5.6) 

n,m>0 

The 3x3 flavour matrix U contains the pseudoscalar meson fields |^ and 

V^U = d^U - tiv^ + a^)U + iU{v^ - a^) . (5.7) 

Here v^, are external vector and axial- vector fields, while the field x = 
2 BoA4 + ... contains the current quark mass matrix Ai. The photon field is 
contained in v^. 

The field strengths R{L)fj_^ are constructed out of the fields r(Z)^ = f^ + (— )a^. 
The covariant derivatives act on R{L)niy as follows 

V^R{L)^, = d^R{L)^, - z[r(/),, R{L)^,]. (5.8) 



The form of the terms in (|5.6|) has been chosen such that there are no vertices 
with three or more photons interacting with pions generated. The first line is the 
lowest order CHPT Lagrangian. The second line contains one- or two-photons 
couplings to pseudoscalar mesons while the last line contains only two-photons 
couplings to pseudoscalar mesons. 

The vertex for a charged pion with incoming momentum p and a photon with 
outgoing momentum k and polarization is given by 



|e| (2p - kf \ g^i, + (k'^g^u - k^k^) ^ E (~^^) 



(5.9) 



The first term is the lowest order vertex. With the choice 

2 \M^J 



(5.10) 



this reproduces the phenomenological complete VMD behaviour for the 'y*P^P 
vertex. But notice that any VMD-like M^{k^) / {M^{k'^) - k'^) behaviour, with 
My(0) 7^ 0, as the one obtained in the ENJL model can be reproduced by choos- 
ing the a„ appropriately. I.e. 

We now turn to the 7*7*P+P^ vertex. If the two photons outgoing momenta 
are ki and /c2 and their polarizations and respectively, that vertex is given 
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by 



2ie g^y 

J IT 



n>0 



J IT 



n>0 

-2z^ (^^fi-^" - ki^k^^ [klgau - k2ak2u 



X 



n,m>0 



Here we see that, keeping gauge and chiral invariance fully satisfied, this vertex 
is rather unconstrained. The choice bnm = reproduces the HGS model used in 



IS, Bl^] with a = 2. The difference between that model and the complete VMD 



model only starts at order p in the chiral counting. 

A large number of other choices are however possible. In particular the choice 

2 

(5.13) 

/ TT 

reproduces the complete VMD amplitude for 7*7*P+P~ used in this work. Again 
depending on the coefficients a^, one can have the ENJL vector meson propagator 
or any other one like the phenomenological complete VMD mentioned above. For 
comparison we use both, see Section |^. It is also possible to add chiral invariant 
terms that will produce a direct dependence on ki ■ k2- This last possibility is 
realized by adding terms like 

(V'^'^V'^V^'L^^) (v^'^V,3'D''Rau) U. (5.14) 

The numerical results for the pseudoscalar mesons loops contribution are dis- 
cussed in Section 0. 



6 Numerical Results 

In this section we give the numerical results for the low energy calculation of the 
hadronic light-by-light contributions to presented in Sections | and ^. 

We first analyze the result for the large Nc limit calculation including the 
effects of the U(1)a anomaly as explained in Section ^. Since we are dealing with 
a low-energy model, as mentioned before, it is necessary to study the dependence 
on a high-energy cut-off /i on the vector legs' momenta. 

For the seven dimensional integral of the pure four-point function (or one 
constituent quark loop) contribution, we used a statistics of 20 iterations with 
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Cut-off 


n,, X 1 O^'' from 


n,. X 1 O^'' from 


n,. X 1 0^'^ from 


a, X 10^° 


u 


Constituent 


Pseudoscalar 


TT*^, rj and rj' 


(GeV) 


Quark 


Exchange O^N^) 


Exchanges 


Sum 




in Figure 3(a) 


in Figure 3(b) 


O(iVe) + U(1)a 




0.5 


0.78(0.01) 


-14.2(0.1) 


-4.8(0.1) 


-4.0 


0.7 


1.14(0.02) 


-19.4(0.1) 


-6.8(0.1) 


-5.7 


1.0 


1.44(0.03) 


-24.2(0.2) 


-9.0(0.1) 


-7.6 


2.0 


1.78(0.04) 


-33.0(0.2) 


-12.6(0.2) 


-10.8 


4.0 


1.98(0.05) 


-39.6(0.6) 


-15.0(0.2) 


-13.0 


8.0 


2.00(0.08) 


-46.3(1.5) 


-17.6(0.4) 


-15.6 



Table 2: Results for the order Nc constituent quark loop and pseudoscalar ex- 
change hadronic light-by-light contributions to in the ENJL model. 



10^ points in the Monte Carlo routine VEGAS, while for the two muon loops 
five dimensional integral of the three-point-like function contributions we used a 
statistics of 20 iterations with 5000 points in the same Monte Carlo routine. This 
statistics is equivalent to the one used in the seven dimensional integral case. For 
the two- and three-point functions needed in these three-point-like contributions 
we used Gaussian integration with an accuracy of 10~^. 



6.1 Pure Four-Point Function 

In Table 1^ we have listed the leading hadronic light-by-light 0{N(^ contributions 
to a^, i.e the pure four-point function in the second column and the pseudoscalar 
exchange three-point-like function in the third column, as a function of the cut- 
off together with the errors quoted by VEGAS. Since the integrand is rather 
irregular, this error estimate is somewhat on the small side (see also [0) and 
will be largely superseded by the error in our final result. 

For the bare constituent quark loop, the result only stabilizes at a rather 
high value of /x. For instance, for a bare quark loop with a constituent quark 
mass of 300 MeV, the change between a cut-off of 2 GeV to a cut-off of 4 GeV 
is still typically 20%. The change from 0.7 GeV to 2 GeV is typically a factor 
of 1.8. These are the results quoted in column 4 of Table [4.1| . The changes 
for our more realistic ENJL model four-point function can be judged from the 
results in Table 0, column 2. It still has a significant change between 1 and 
2 GeV cut-off. This invalidates the use of any low energy model to calculate 
accurately the complete hadronic light-by-light contribution to a^. The bulk of 
these contributions does not come from the dynamics at scales around the muon 
mass as it is often stated. This also explains the rather high sensitivity to the 
damping provided by the vector two-point functions as seen in |jl6[ . Mostly due to 
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Cut-off 
/i 

(GeV) 


X 10^° 
ENJL 


X 10^° 
Point-Like-VMD (tt") 


X 10^° 
V 


X 10^° 

7]' 


0.4 


-2.84(2) 


-2.70(1) 


-0.425(1) 


-0.266(1) 


0.5 


-3.70(3) 


-3.46(2) 


-0.616(2) 


-0.399(2) 


0.7 


-5.04(4) 


-4.49(3) 


-0.923(3) 


-0.631(2) 


1.0 


-6.44(7) 


-5.18(3) 


-1.180(4) 


-0.847(3) 


2.0 


-8.83(17) 


-5.62(5) 


-1.37(1) 


-1.03(1) 


4.0 


-10.51(37) 


-5.58(5) 


-1.38(1) 


-1.04(1) 



Table 3: The vr^ exchange contribution to for the ENJL and the point-hke 
Wess-Zumino vertex, damped with two vector propagators for the 7r°, t] and rj'. 



its electric charge and heavier mass, the contribution of the strange quark flavour 
is much smaller (around 0.04- 10~^°) than that of the up and down quarks shown 
in the Table ^ This value is within the quoted VEGAS error for up and down 
quark contributions. In Section ^ we give an estimate of the intermediate and 
high energy one constituent light quark loop contributions and the heavier quark 
flavours contributions. 



6.2 Pseudoscalar Three-Point-Like Function 

For the three-point-like function contributions we have done the same study of the 
cut-off dependence as for the four-point function contribution. In particular we 
find that at large N^. the contribution of the pseudoscalar exchange is more than 
one order of magnitude larger than the others. The reason this contribution is so 
different can be traced back both to the presence of two flavour anomaly vertices 
and the CHPT counting. It therefore deserves more attention. In fact, the 
pseudoscalar exchange has important next-to-leading corrections from the effects 
of the U(l)^ anomaly that leave the 7r° exchange as the dominant contribution 
to a|^gi^t-t'y-i'si^t. First we give in column 3 of Table ^ the result strictly to leading 
order in 1/Nc from the u, d and s flavours. 

We have taken into account the effects of the U(l)yi anomaly by using the 
physical 7r°, rj and rj' mass eigenstates as propagating states instead of the mm, dd 
and ss in the large N^. limit. This includes the main effect of the U(l)^ anomaly 
which is in the differences of the masses of the pseudoscalar rj and rj' mesons. 
In the 7r°, rj and rj' basis, the contribution from the pion intermediate state has 

r n 2 

a charge factor ((2/3)^ — (—1/3)^) / \/2 compared to a single quark of charge 
one. We thus multiply the ENJL pseudoscalar exchange result of column 3 in 
(^ by this factor to obtain the vr'' contribution from the ENJL model listed in 
column 2 of Table ^ 
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The 7] and t]' contributions we cannot estimate directly within the ENJL 
model. What we have done is the following. Since the main effect is in the 
propagator of the exchanged pseudoscalar meson, the ratio of the tt" exchange 
contribution to the 7] or 77' contribution has to be in good approximation model 
independent. So we have taken the ratio from the point-like Wess-Zumino La- 
grangian with full vector meson dominance and multiplied the ENJL tc^ contri- 
bution by them to get the 1] and 1]' contributions. So to get the fourth column in 
Table ^ we sum the last three columns in Table multiply by the 2nd column 
and divide by the third one. This is our estimate for the combined 7r°, 77 and t]' 
contributions. For the calculation of the last three columns in Table ^ we have 
used P°77 couplings such that the experimental decay rates P° 77 are repro- 
duced and a vector meson mass of 0.78 GeV. That gives ratios that vary from 
16% at /i = 0.4 GeV to 25% at /i = 4.0 GeV for the ratio of the r] contribution 
to the 7r° one and from 10% at /i = 0.4 GeV to 19% at /z = 4.0 GeV for the ratio 
of the 77' contribution to the 7r° one. 

We find for the pseudoscalar result less stability at high values of the cut- 
off fi than for the quark- loop contribution. Although the change from 0.7 GeV 
to 2 GeV is also around 1.8, the stability is worse for cut-off values above 4 
GeV. Notice also that the error from the integration routine VEGAS is larger 
for these values of the cut-off. The poor stability in the pseudoscalar exchange is 
mainly due to the subtraction terms we need to obtain the correct SU(3) flavour 
anomaly. We shall give in Section |^ an estimate of the intermediate and high 
energy contributions to the pseudoscalar exchange term. We finally give in the 
fifth column of (0) the sum of the second and fourth columns. 

6.3 Other Three-Point-Like Functions 

Both scalar and axial-vector exchanges in three-point-like function contributions 
are much smaller than our final error. Their results for up, down and strange 
quark flavours are in Table ^ The scalar contribution has obviously stabilized. 
The axial-vector one has large cancellations and becomes numerically unstable 
for a cut-off of 8 GeV. We have therefore not quoted the values for this cut-off. 

6.4 Pion and Kaon Loops {0{1) in l/N^ Contributions) 

The results for the dominant contributions of order 1 in 1/Nc are in Table |^. We 
have saturated this contribution by the physics of pion, kaon and rho mesons as 
explained in Section ^ Therefore, we need to verify if these contributions really 
saturate for energies below the axial-vector mass. For this, we have studied the 
cut-off dependence by varying the Euclidean cut-off /i. For these contributions 
we used a statistics of 10 iterations with 10^ points in the Monte Carlo routine 
VEGAS. As can be seen in Table ^, the charged pion loop contribution saturates 
around 2 GeV, while the kaon loop contribution saturates around 4 GeV. From 
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Cut-off 


X 10^° from 


X 10^° from 


X 10^° 




Scalar 


Axial- Vector 


(GeV) 


Exchange 0{Nc) 


Exchange 0{Nc) 


Sum 




in Figure 3(b) 


in Figure 3(b) 




0.5 


-0.22(0.01) 


-0.05(0.01) 


-0.27 


0.7 


-0.46(0.01) 


-0.07(0.01) 


-0.53 


1.0 


-0.60(0.01) 


-0.13(0.01) 


-0.73 


2.0 


-0.68(0.01) 


-0.24(0.02) 


-0.92 


4.0 


-0.68(0.01) 


-0.59(0.07) 


-1.27 



Table 4: Results for the order Nc scalar and axial- vector exchange hadronic light- 
by-hght contributions to in the ENJL model. 



Cut-off 


X 10^° from 


a^t X 10^° from 


X 10^° 


fi 


Pion Loop 


Kaon Loop 


(GeV) 


in Figure 3(a) 


in Figure 3(a) 


Sum 


0.5 


-1.20(0.03) 


-0.020(0.001) 


-1.22 


0.6 


-1.42(0.03) 


-0.026(0.001) 


-1.45 


0.7 


-1.56(0.03) 


-0.034(0.001) 


-1.59 


0.8 


-1.67(0.04) 


-0.042(0.001) 


-1.71 


1.0 


-1.81(0.05) 


-0.048(0.002) 


-1.86 


2.0 


-2.16(0.06) 


-0.087(0.005) 


-2.25 


4.0 


-2.18(0.07) 


-0.099(0.005) 


-2.28 



Table 5: Results for the order 1 in the 1/Nc expansion charged pion and kaon 
loops hadronic light-by- light contributions to a^. 



27 



Cut-off 




fi 


X 10^° 


(GeV) 




0.5 


-5.5(0.1) 


0.7 


-7.8(0.1) 


1.0 


-10.2(0.1) 


2.0 


-14.0(0.2) 


4.0 


-16.6(0.2) 



Table 6: The contribution to from the low energy domain. The error is the 
one of the various contributions given by VEGAS and added in quadrature. 



Table ^, we see that the change between the result at 1 GeV and the result where 
it stabilizes is less than 20% so we conclude that the approximation we are doing 
works to this accuracy which is good enough in view of our final uncertainty, see 
Section ^. The intermediate and higher energy contributions for this case are 
discussed in the next section. The results in Table ^ are obtained using ENJL 
vector mesons propagators for the vector legs. We have also used vector meson 
propagators with a constant vector mass of 768 MeV and in this case the charged 
pion plus kaon loop contributions to a^^^^^-^y^^^s^^ saturate earlier (at fi = 0.8 
GeV ) with a value around —1.65 ■ 10~^°. 



The HGS model with a=2 was used in Ref. [|I^, ^ to calculate this contribu- 
tion. We see from ( |3.4| ) that the HGS in the non-anomalous sector and for a = 2 
has a wrong high energy behaviour when matching QCD in the 7r+ — vr" mass 
difference. See the negative u"^ correction to the logarithmic behaviour there. 
This also tends to lower the contribution to too much when vector mesons are 
added. This is important in the charged pion and kaon loop contributions where 
the P"''P~7*7* has an unknown high energy behaviour. We have adopted the 
criterion of using a complete VMD model inspired by the 0{Nc) ENJL model. 
As shown in Section ^, this does not break any Ward identity. This choice has, at 
least, a good high energy behaviour for two-point functions, e.g. Weinberg Sum 



Rules are fulfilled. This is not true for the HGS model with a = 2. See also 34 



where the tt^ — n mass difference is calculated within this model. The result of 
the HGS model can however not be excluded with these arguments. 



6.5 Sum of Low-Energy Contributions 

Adding the contributions calculated before we get the low energy contribution to 
as a function of the cut-off fi. These are the final results for the low-energy 
contribution estimated within the simplest version of the ENJL model. We will 
present other estimates in the next section. The results are in Table ^. The ENJL 
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model we have used is a good low energy hadronic model which works within 20% 
up to energies (0.4~0.6) GeV depending on the channel. We observe from the 
results in Table |^ that higher energy contributions are certainly not negligible. 
The estimation of those contributions is the subject of the next section. We 
conclude from this section 

^light-by-iight^^ = 0.5 GeV) = -5.5(1.5) ■ 10-^°. (6.1) 

The error includes five times the integration error from VEGAS and the estimated 
model dependence added in quadrature. This is rather small for the dominant 
pseudoscalar exchange contributions, as is shown by the small changes in the 
various models presented in the next section. At these energies the main error is 
from the model dependence of the pseudoscalar meson 0{1) contribution which 
we estimate to be 0.8 ■ 10^^° to cover the results in [|1^, |21|] as mentioned above. 



7 Intermediate and High Energy Contributions 

In this section we estimate the hadronic contributions from intermediate and high 
energy regions to a|^s'^*~^^~^'s'^* . Here we are already outside the applicability of 
the ENJL model that we have only used for the low energy region. 

We want to make a general comment regarding the use of the ENJL factor 
Myl—p^)/ {My{—p'^) — p^) at large Euclidean scales. If we naivelyQ send the Eu- 
clidean = —p^ — > cxo, then this factor becomes 1. Therefore the photon x vector 



propagator in Eq. (O) goes to zero as l/Q^ in the ENJL model while in a VMD 
model it goes to zero as My/Q^. This difference will not affect very much the 
calculation since at low energies where we apply the ENJL model, both vector 
meson propagators behave very similarly. 



7.1 Pure Four-point Function 

In the case of the constituent quark loop contribution one can still obtain an 
estimate of the higher energy contributions, e.g. by mimicking the high energy 
behaviour of QCD by a bare constituent quark loop with a mass of about 1.5 GeV. 
This gives only an additional correction of 0.24(1) ■ 10~^°. Using the results of [16 
this scales like 1/M^ with M the quark mass. We have checked this behaviour 
as well. Here the mass of the heavy quark acts as an infrared cut-off so that this 
heavy bare quark loop is mimicking the QCD behaviour for a massless quark with 
an IR cut-off around 1.5 GeV. We take this number both as the value and the 
uncertainty due to the high energy region contribution for the three light flavours. 
In fact we can simply assume that n = M and add the ENJL contribution up to 
the scale to the bare quark-loop contribution with mass M in the loop. This 



^Notice that this cannot be done in the ENJL model since is only valid for energies \p\ << A. 
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leads to a total light-quark contribution of 2.2, 2.0, 1.9 and 2.0 (times 10^^°) for 
a matching scale /i = M of 0.7, 1, 2 and 4 GeV, respectively. 

We estimate the charm quark contribution with a bare quark loop. If we damp 
it with cc meson dominance propagators in the photon legs it will be somewhat 
smaller. This contribution is also small (about 0.2 ■ 10~^°). Therefore we obtain 
for the total contribution from the pure four-point function part: 

a^(Quark - Loop) = 2.1(0.3) ■ 10"^°. (7.1) 



7.2 Pseudoscalar Three-Point-Like Function 

Here we will only discuss the estimate of intermediate and high energy con- 
tributions from the pseudoscalar exchange since this is the dominant hadronic 
contribution to al:g^*-^y-"g^*. As seen in Section 6 all others are much smaller 
and changes there will not affect our result significantly. 

This contribution can be seen as the convolution of two P^'j*'-f* vertices with 
both and photons off-shell and in the Euclidean region (see Eq. (|4.10|) ). Let 
us summarize what we know about this form factor. The 7r°7*7 form factor has 
been measured at CELLO and CLEO-II |^3| for values of the Euclidean in- 



variant photon mass above (0.8)^ GeV^ and below (2.8)^ GeV^. These are 
the data points in Figure As we mentioned before, we know that in QCD the 
7r°7*7* vertex has a 1/Q^ asymptotic behaviour when one of the photon legs Eu- 
clidean momentum is very large |^ . This is supported by the phenomenological 



analysis of J/\l/ decays in the same reference, where one finds that the asymptotic 
behaviour predicted by QCD works reasonably well from scales around the J/\l/ 
mass. 

The chiral anomaly only fixes the P°77 form factor at 0(p'^), and this is ful- 
filled by the ENJL model form factor. To estimate the pseudoscalar-exchange 
intermediate and high energy contribution to a^, what we have done is to find 
a phenomenological parametrization that interpolates between the ENJL form 
factor, which is supposed to work well below 0.5 GeV, and the measured 7r°7*7 
form factor for Euclidean energies above 0.5 GeV and with its asymptotic be- 
haviour predicted by QCD at large Euclidean momentum. Notice that what we 
really need is the 7r°7*7* form factor as we said above. Unfortunately no data 
are available for this form factor. 

Different parametrizations give quite different contributions to for cut-offs 
larger than 0.5 GeV. A lower limit would be the results from the point-like Wess- 
Zumino vertex without vector meson dominance. This assumes that the pion 
remains pointlike at all relevant scales. Its contribution to is given in column 
2 of Table |^ for the vr" exchange. The logarithmic behaviour as a function of the 
cut-off is clearly visible here. The results from the ENJL model we used for the 
low energy contribution are in column 2 of Table |^. As can be seen even at a 
cut-off of 0.5 GeV the effect of the damping in the ENJL model is already very 
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Cut-off 


X 10^° 


X 10^° 


X 10^° 


X 10^° 


X 10^° 


/i 






Point-Like- 


Transverse- 


Transverse- 


(GeV) 


Point-like 


ENJL-VMD 


VMD 


VMD 


VMD 


0.5 


- 4.92(2) 


- 3.29(2) 


- 3.46(2) 


-3.60(3) 


-3.53(2) 


0.7 


-7.68(4) 


- 4.24(4) 


- 4.49(3) 


-4.73(4) 


-4.57(4) 


1.0 


- 11.15(7) 


-4.90(5) 


- 5.18(3) 


-5.61(6) 


-5.29(5) 


2.0 


-21.3(2) 


- 5.63(8) 


- 5.62(5) 


-6.39(9) 


-5.89(8) 


4.0 


-32.7(5) 


- 6.22(17) 


- 5.58(5) 


-6.59(16) 


-6.02(10) 



Table 7: The 7r° exchange contributions above n = 0.5 GeV in various 
parametrizations of the 7r°7*7* vertex that fit the data for the 7r°77* vertex. 



important. 

Both these parametrizations (ENJL and point-like Wess-Zumino vertex) do 
not fit the measured data points for the 7r°7*7 form factor above (0.5~0.6)^ GeV^ 
for the Euclidean photon invariant mass, see Figure ^. 

The simple point-like Wess-Zumino vertex plus VMD fits the 77^77* data at 
high off-shellness reasonably well (see Figure This we call Point-Like- VMD 
parametrization and the results of the vr'' exchange contribution to are quoted 
in column 4 of Table 0. This parametrization is suppressed in both photon prop- 
agators and could be too suppressed for both photons far off-shell. It can thus 
be taken as an upper limit for the 7r° exchange of this type of contribution for 
high Euclidean momenta. As another prescription to interpolate between the 
low-energy ENJL form factor and the measured form factor we can multiply 
the subtraction terms in the ENJL model P°7*7* form factor with ENJL vector 
meson propagators, My{—p^)/ {My^—p"^) —p"^), in all photon legs, this we call 
ENJL-VMD parametrization. This is somewhat too hard at very high energies 
but reproduces the form factor data at intermediate momenta reasonably well. 
The tt" exchange contribution to for this parametrization is in column 3 of 
Table 0. We observe that both Point-Like- VMD and ENJL-VMD parametriza- 
tions give numerically very similar contributions to a^. Another parametrization 
of the 7r°7*7* amplitude that also fits the 77^7*7 form factor data is the following 

J"''iPuP2) = ^^^e^'^"W2/.^(p?,P2, q') (7.2) 

with 

M2 



HpIpIq') = {l-9A{~q' 



M2-(pi-p2)2 



31 



0.5 



1.5 
mu (GeV) 



Figure 4: The ir^'j*'^ form factor. At high energies the curves from top to bot- 
tom are: ENJL, ENJL-VMD, Point-hke-VMD and Transverse-VMD. The data 
points are obtained from scahng the TPC/Two-Gamma rj resuhs |0 (circle), the 
CELLO results (crosses) and CLEO-II (squares). The horizontal error 
bars are 0.1 GeV for the TPC/Two-Gamma and CELLO data and 0.05 GeV for 
the CLEO ones. 
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mU-pD - pI m^-pD - pr ^^-^^ 

The form factor F{pl,p\,ql) is defined in Appendix 0, Eq. ( p.6|) , gA{,—'f') and 
M^(-p2) in Ref. |8|. It fits the data when the " mass" M varies within 0.6 < 
M < 1.4 GeV.The results for this parametrization that we call Transverse- VMD 
are in columns 5 and 6 of Table for a "mass" M of a/2 and 1 GeV, respectively. 

In general, we obtain that the parametrizations that fit the data also match 
the ENJL low energy result at 0.5 GeV reasonably well (compare Tables |] and 
1^). This supports the low model dependence error for the lower than /x = 0.5 
GeV energy domain pseudoscalar exchange contributions. 

We estimated the effects of rj and r]' in the same way as was done in the 
previous section. 

Let us analyze the pseudoscalar exchange contributions from scales higher 
than 4 GeV. The QCD behaviour predicted in |^ goes like l/Q^, where 



is the invariant Euclidean mass of the off-shell photon. This dependence 
suppresses the high energy contributions more than the point-like Wess-Zumino 
vertex damped by complete VMD propagators. In this last case we see that the 
contributions from energies above 2 GeV are negligible. So we consider those 
from above 4 GeV negligible. 

We do not take the simple ENJL model as given in the previous section but 
the parametrizations that fit the form factor data. Since their results are very 
similar, we average the result for the pointlike with VMD factors, the ENJL- VMD 
and the transverse- VMD with M = 1 GeV at 4 GeV. This is our final result for 
the pseudoscalar exchange: 

(Pseudoscalar - Exchange) = -8.5(1.3) ■ 10~^°. (7.4) 

Here the error is estimated as about 15%. This includes all the models mentioned 
except the pure ENJL result. 



7.3 Other Three-Point-Like Functions 

Since the contribution of the scalar exchange is smaller than the final error of the 
dominant pseudoscalar meson exchange, we quote for it 

a/,(Scalar) = -0.68(0.2) ■ 10"^°. (7.5) 

without making any futher estimate of the suppressed higher energy contributions 
in this case. There are in any case no pointlike subtractions here that would have 
produced large possible changes. 

For the axial-vector exchange contribution we have done an analysis similar 
to the pseudoscalar one. We analyzed the simplest ENJL model in Section ^ 
and we can do the analog of the ENJL- VMD model as well. This corresponds to 
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Cut-off 


a,, X 10^° from 


/i 


Axial- Vector 


(GeV) 


Exchange 0{Nc) 




in Figure 3(b) 


0.5 


-0.04(0.01) 


0.7 


-0.06(0.01) 


1.0 


-0.10(0.01) 


2.0 


-0.15(0.01) 


4.0 


-0.35(0.04) 



Table 8: Results for the axial- vector exchange in the ENJL-VMD parametriza- 
tion. 

removing the VMD-factors in Eq. ( [B.15| ). The results for the ENJL-VMD form 
factor is in Table (^. Taking into account that there are cancellations which 
cause VEGAS to underestimate the error at the scale of 4 GeV, we take as the 
final result for the axial exchange 

a^(Axial - Vector) = -0.25(0.1) ■ 10"^°. (7.6) 

7.4 Pion and Kaon Loops {0{1) in 1/Nc Contributions) 

As explained in Section]^, we have saturated the C(l) contributions with charged 
pion and kaon loops modulated by vector meson propagators for the vector legs. 
Since our model is only valid for energies below the axial-vector mass, we take 
the difference between the result at /i = 1 GeV and where it stabilizes as an 
estimate of contributions from resonances heavier than 1 GeV running in the 
loop. These contributions are higher order in the chiral counting and suppressed 
by inverse powers of the mass of these resonances. Therefore, we take them as an 
estimate of the intermediate and high energy contributions for this 0{1) in the 
1/A^c expansion contribution. For the 0{1) in the l/N^ expansion contributions 
to we thus quote 

l^^ught-by-iight^^^^^ = -1.9(1.3) ■ 10"^° , (7.7) 

where we have taken as central value the result at = 1 GeV and as error the 
high energy contribution as estimated above plus five times the VEGAS error 
added in quadrature. There is an extra 0.8 added linearly to the error because of 
model dependence, see comments in Section ^ about its origin. It also includes an 
educated guess of the 0(1/A^c) corrections and the rest of the 0{1) contributions. 
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8 Discussion of Results and Conclusions 



Our final estimate for the hadronic light-by-light contributions to and main 
result of this work is the sum of the partial contributions in Eqs. ([7.1|), ( |7.4| ), 
(|7l) and (O), 



^light-by-iight _ _9_2(3.2) . 10-1° . (8.1) 



The error is obtained by adding linearly the error of each contribution. This 
is because we are essentially using the same model for all contributions so the 
error is likely to be in the same direction for all contributions. This results in a 
35% error which we believe takes adequately into account the model dependence 
error of this calculation. This result improves and substitutes the one in Ref. 



There we took as first estimate the result in Ref. [19| for the 0{1) in the 
1/Nc expansion contributions. Here, we have given an estimate for it and a more 
detailed analysis of the high energy contributions has been performed. We also 
performed a more detailed study of the rj and t]' effects. 

We want now to present the result in ( ^.1| ) by explicitly splitting the different 
contributions. First, separating the high and the low energy contributions we get 

^light-by-iight _ Q _ 4 . 10-10, (8.2) 

where the first number is the lower than n = 0.5 GeV contribution (for definite- 
ness in the ENJL-VMD model, the others range from —5.0 ■ 10~i° to —5.5 ■ 10~i°) 
and the second the higher. Here one can see that the intermediate and higher en- 
ergy contributions are certainly not negligible and a^^s^^~^^~^^s^^ does not saturate 
at low energies as assumed in |jT5|, [16 1. 



It is also interesting to see how well the l/N^ expansion works in this case. 
The leading l/Nc result is about — 22- 10~i°, so the 1/Nc correction is about 50%. 
Notice, however, that most of it is from the U(l)^ anomaly contribution which 
does not appear at order The 1/Nc expansion works thus OK. 

Finally, we see that all contributions except for the pseudoscalar one cancel 
to a large extent. The part from the pseudoscalar alone is —8.5 ■ 10~i°. We see 
clearly the dominance of the pseudoscalar exchange. Notice however the large 
cancellations occurring between the four-point- like function contribution and the 
pion and kaon loops ones, namely (2.1 — 1.9) ■ 10~i°. This is possible since the 
pion loop is suppressed by 1/Nc but dominant in the chiral counting while the 
other is suppressed by the chiral counting but leading in 1/Nc. 



we 



Since the works and [|T^] have been updated and/or corrected by 
refer to this last one for the comparison with other works. The authors of [jl9 
get as final result 

light-by-light _ _5.2(1.8) ■ 10-1°. (8.3) 



a 



The main differences from the results in [jT9[ are the following. The 77', scalar and 
axial-vector exchange contributions were not included there, this amounts to 

(-1.11 - 0.68 - 0.25) ■ 10-1° = -2.04 ■ 10-i°. (8.4) 
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They have to be certainly included. 

Another difference comes from the different estimation of the pseudoscalar 
meson loop 0{1) in 1/Nc contribution. We have essentially used complete VMD 
for the P~^P~'j*'j* vertices for different reasons, see Sections ^ and ^ for a detailed 



explanation. The numerical difference with [11^] for this contribution amounts to 



-1.45-10-1". (8.5) 

We have taken into account this model dependence by adding linearly an extra 
factor to the error in ( |8.1| ). 

The rest of the numerical discrepancy is small (—0.5 ■ lO'^'' ) and due to 
several causes: simplifications in [jl9l, VEGAS numerical uncertainty, .... Its 



smallness is gratifying and reflects the low model dependence of the rest of the 
contributions. 

Our calculation establishes that the contribution to from light-by-light 
scattering is negative and relatively large. It is one half of the electroweak cor- 
rections 0]. This result is between two and three times the aimed experimental 
uncertainty at BNL. Although we believe our error estimate is conservative, it 
has an unsatisfactory uncertainty that will be difficult to reduce because of model 
dependence. This is mainly in the pseudoscalar exchange and the pseudoscalar 
meson loop contributions. Despite this uncertainty, the estimate in ( fj.lD is still an 
important theoretical result for the interpretation of the muon g — 2 measurement 
at the planned BNL experiment. 

Adding the theoretical calculations of the Standard Model contributions to 
in Eqs. (|r|), (|OD, (|r|), (g^), and gives the following present 



theoretical estimate for the muon g — 2, 

a^^ = 11659182(16) ■ 10-1°, (8.6) 

where the quoted errors for different contributions are added in quadrature. Using 
for the full photon vacuum polarization insertion in the electromagnetic muon 
vertex the result in ( |1.6|) instead of the one in ( |1.5|) p|, gives 



a*'' = 11659168(11) ■ 10-1°. (8.7) 
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A Construction of W""^ 

In this appendix we give the general Lorentz structure of n'"^"'^(pi,p2,P3) defined 
in Eq. ( p.3| ). See Fig. |1] for definition of the momenta. This four-point function 
can be decomposed by using Lorentz covariance as follows 

+ 0=^(^1,^2,^3)^/^^^?^" 

+ Il'''{pi,P2,P3)9'^p',pt + ^''HPl,P2,P3)9''yyk 
+ n'^HpuP2,P3)9''y,Pt + n^''(Pl,P2,P3)^"MPfc 

+ ff^''-™(pi,P2,P3)pfp>&:;, (A.i) 

where fc,m = 1, 2 or 3 and repeated indices are summed. There are in total 
138 Ll-functions. Not all of them are independent since they are related by Ward 
identities. In fact, when all possible permutations of the vector legs contributing 
to W^°''^{pi,p2,P3) are summed, U(l) gauge covariance leads to the following 
Ward identities 

PluIl"""'{Pl,P2,P3) = P2aIl"""'{Pl,P2,P3) = 
P3/3n^^"^(pi,P2,P3) = qpUP^"^^{p^,p2,P3) = . (A.2) 

The extensive use of these relations makes possible to express the amplitudes 
n^'^'^(pi,P25P3) and Il^^^' "'^^^{pi,p2,P3) in terms of Il''^''"^{pi,p2,P3). The fact 
that we can write down everything in terms of ir^^'^{pi,p2,P3) is again telling 
us that the light-by-light scattering contribution to finite quantity. We 

just have to calculate the UV safe W^^"^{pi,p2,P3) amplitudes using, for instance, 
a cut-off regularization scheme like proper-time. This regulator introduces the 
physical cut-off A of the ENJL model, see Section |^. 

The amplitudes W^^''^{pi,p2,P3) are not the minimal set of independent am- 
plitudes. We can still reduce it further with some additional Ward identities. 
However, we shall not use all of them and leave some ones as checks on the 
resulting W^°'^{pi,p2,P3). Since the quantity we need to compute is the anti- 
symmetric part of M'^^(O) in Eqs. ( |2.6|) and ( p.7|) , which contains the derivative 
Qf n^^"'^(p-,^^p2,P3) with respect to p3x at ps = 0, we can reduce the number of 
needed amplitudes to n^-''^'"(pi,p2,P3), n*^^"'(pi,p2,P3), ^^^^^iPi,P2,P3) and the 
derivatives oiW^''^{pi,p2,P3) — Il''^^'^{pi,p2,P3) with respect to p3\ at p3 = 0. Here 
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i,j,k,m = 1 or 2, so that we need 32 functions. This is the set of amphtudes 
which we will use in all our calculations. 



B Three-Point Functions 

In this appendix we give the three-point functions needed in Section ^. Barred 
three-point functions are the one-constituent-quark loop ones. From Eqs. ([4.7|) , 
( ^.10| ) and ( |4.11| ) we see that we only need barred three-point functions, therefore 



we only give the explicit expression for them. The full three-point functions can 
be obtained using the methods explained in pS| in a straightforward manner. 



We only will give the contribution to the three-point function given by the clock- 
wise orientation of the internal quark lines. The other orientation is taken into 
account in the permutation of the external vector legs, see Section ^. 
Let us start with the SVV three-point function. This is defined as 

nfrbi,P2) = ^'J d'x J d'ye'^^-^^^-y^ (0|T (s^^(o)v^'(x)v;™(y)) |0), 

(B.l) 

with V^^{x) = [qi{x)'y^qj{x)] and S''''{y) = —[qk{x)q{x)i]. Summation over colour 
between brackets is understood and latin indices are flavour indices. Owing to 
Lorentz covariance this three-point function can be decomposed as follows 

+ U3{pl,p2)p2^,Plu + 'n^4{pl,P2)p2f,P2u + 'n^5{Pl,P2) gflu ■ 

(B.2) 

We shall use the Ward identities for three-point functions to write down Il^{pi,p2) 
in terms of the W{pi,p2), i = 1, ■ ■ ■ ,4 amplitudes which are UV finite. We com- 
pute the corresponding barred functions Il\pi, P2) , i = I, ■ ■ ■ , 4: with the standard 
Feynman parametrization technique and using proper-time regularization. This 
regulator introduces the physical cut-off A, see Section ^ The VVS three-point 
function can be obtained from this one using the identity 

K^'{PuP2) = njr(-(Pi +P2),Pl) . (B.3) 
The anomalous PVV three-point function is defined as 

Il^nPi,P2) ^ t'Jd^xJd^ye'^^-^+P^-yH0\T{p'^{0)V;:^ix)Vr iy)) |0) 



(B.4) 



with P'^^{x) = [qi{x)i'y5qj{x)]. It was calculated in the ENJL model in Ref. 
With the same notation as there, we get 

-/Vj-; Q Q / 9 9 9 \ ^ 
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with Mi is the flavour i constituent quark mass. Function ( |B.5| ) is the one- 
loop constituent quark (barred function) contribution with clock- wise orienta- 
tion of the internal quark line to the three-point function ( [B.4| ). The function 
F{plplq^) is 



Fiplplq') = l-h{M^Q,0,0,0) + UM^,plplq''), , 



(B.6) 



where Mq is the constituent quark mass in the chiral limit and 
where 



M^{x, y) = M^- x{l - x)pI - y{l - y)pl - 2xyp, ■ p^ , (B.8) 
A is the physical cut-off, see Section ^ and 

^l(e)=e-^ (B.9) 



Notice that in Eq. (p.6| ), we have given the constituent quark mass dependence 
that was not explicit in Ref. p8|. In Eq. we have corrected an obvious 

misprint in Ref. [p^] . 

Again the VVP three-point function can be obtained from the PVV one by 
symmetry 

nJT(Pi,P2) = njr(-(pi + P2),Pi) . (B.io) 

The anomalous AVV three-point function is defined as 

n^r/(pi,P2) = d'x I d'ye^^^-^^^-y^ (o|T (4(o)K''(x)vr"(y)) |0), 

(B.ll) 

with ^^(x) = [qi{x)'y^'y5qj{x)]. The most general Lorentz decomposition in four 
dimensions of this three-point function is 

^7vviPuP2) = te^'"'^{p^pU^{puP2)+P2(sIl2{PuP2)) 

+ te>^''^'p,,P2S Kn3(pi,P2) +P^n4(pi,P2)) 

+ ze''''^'p,,p2s{p';MPi,P2)+P2MPuP2)) . (B.12) 

We have used Schouten identities to eliminate redundant terms. These identities 



can be also used to relate this basis to the one used in . As for the other three- 
point functions not all the amplitudes Ili{pi,p2), i = l,---,6, are independent 
since they are related through Ward identities. In this case there are three Ward 
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identities that reduce the six amphtudes to three independent ones. For the 
diagonal flavour case, the Ward identities for the two vector legs, give 



Il2{pi,P2) = Pin5(pi,p2) +Pl ■P2^6(P1,P2)■ 



(B.13) 



The Ward identity on the axial leg relates the amplitudes 113(^1,^2), 1^4(^1,^2), 
115(^1,^2) and IIq{pi,P2) to the form factor of the Uj^^^ {pi,p2) anomalous three- 
point function in (|B.6D . The explicit expressions for the corresponding barred 
AVV three-point function amplitudes (pulling out the flavour structure factor 
§im§ni^jk^ are 

n3(pi,P2) 

^4{P1,P2) 
^5{PUP2) 



-Uq{pi,P2) 



167r2 
167r2 ' 



8 / dx 



dy 



xy 



167r2 



dx 

1 



dy 



M^{x,y) 



Jo M"^{x,y) 



t^{M\x,y)/A') 
T,{M\x,y)/A') 



8 / dx 







dy 



x(l - x) ~ 



M^{x,y) 



T,{M\x,y)/A') 



:B.14) 



with M'^{x,y) deflned in Eq. ( |B.8| ) and 

M\x,y) 



T,{M\x,y)/A' 



X 



M^{x,y)-Mf 
M? 



1 



(mU-pD - pi){MU-pl) - pD 



M^{x,y)-Mi 



mU-pI)mU-pI) 

^T,{M\x,y)/A'). 



(B.15) 



The function My{—p^) can be found in [||] and each M^^-pf) in ( [BlsD corre- 
sponds to the vector leg flavour numbers in (|B.4D. Again the VVA three-point 



(B.16) 



function can be obtained from the AVV one by symmetry 



n 



VVA 



{Pl,P2] 



n 



AVVf 



iPi+P2),Pi] 
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